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The emission pattern in the V per p — V par plane of Inter- 
mediate Mass Fragments with Z—3-7 (IMF) has been studied 
in the collision 116 Sn + 93 Nb at 29.5 AMeV as a function of 
the Total Kinetic Energy Loss of the reaction. This pattern 
shows that for peripheral reactions most of IMF's are emitted 
at mid- velocity. Coulomb trajectory calculations demonstrate 
that these IMF's are produced in the early stages of the reac- 
tion and shed light on geometrical details of these emissions, 
suggesting that the IMF's originate both from the neck and 
the surface of the interacting nuclei. 
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An intense emission of intermediate mass fragments 
(IMF) at mid- velocity has been put into evidence 
in collisions of heavy ions in the Fermi 
regime, i.e. at bombarding energies from 30 to 50 AMeV. 
For non-central collisions, where the binary character of 
the reaction is preserved, the mid-velocity particles are 



easily identified in the V pi 



Vp ar plots, the velocity 



components being usually defined with respect to the 
separation axis of the two heavy reaction partners. The 
velocity region between target and projectile presents 
an intensity of particles unexpected on the basis of the 
statistical emission from the hot, fully accelerated main 
fragments. The multiplicity of mid-velocity particles is 
usually obtained by subtracting from the total emission 
the contribution of particles statistically emitted by the 
projectile- and target-like fragments (PLF and TLF), 
see ^M. Presently, the isospin composition of the mid- 
velocity light charged particles (Z < 2, LCP) and IMF's 
is highly debated, as it can be related to the composition 
of the emitting source, the so-called neck (see, e.g., B and 
references therein). Equally debated is the issue of the 
mechanism responsible for these mid-velocity emissions, 
namely whether they are of dynamical origin or due to 
an exotic statistical process, induced by perturbations of 
the Coulomb field §. 

It was shown by the INDRA collaboration Q| that, for 
a given bombarding energy, the largest ratio (up to »3) 
of mid-velocity to statistical IMF's is found in non-central 
collisions. However, due the relatively high thresholds for 
TLF detection of the existing 47r detectors, up to now the 
very peripheral collisions could not be reliably selected 
and the mid-velocity emission of IMF's was never cleanly 
isolated from the accompanying statistical emission, thus 



preventing a detailed study of its characteristics. 

In this Letter we present experimental results which, 
thanks to an efficient selection of peripheral reactions, 
make it possible to clearly isolate the mid-velocity emit- 
ted IMF's from the statistical component. The obtained 
selection allows -for the first time- to identify reactions 
where almost all IMF's concentrate at mid-velocity and 
to carefully study their emission pattern. Preliminary 
results were presented in Ref . [|| . 

The experiment (Florentine Initiative After Supercon- 
ducting Cyclotron Opening, FIASCO) was performed at 
the Laboratori Nazionali del Sud in Catania (Italy) . Tar- 
gets of 93 Nb («200 fig /cm 2 thick) were bombarded with 
a 29.5AMeV pulsed beam of 116 Sn of <0.1 nA intensity 
and <1 ns time resolution. The setup basically consisted 
of 24 position-sensitive Parallel Plate Avalanche Detec- 
tors (PPAD) @0 covering w 70% of the forward solid 
angle. They measured, with very low thresholds (< 0.1 
AMeV), impact time and position (FWHM resolution 
of 600 ps and 3.5 mm, respectively) of heavy reaction 
products (A >20). From the velocity vectors, primary 
(i.e. pre-evaporative) quantities were deduced, event-by- 
event, with an improved version p"2| ] of the Kinematic 
Coincidence Method (KCM). 

With respect to previous experiments |^3|,Q, the 
setup included 160 phoswich scintillators mounted be- 
hind most of the PPAD's. They detected LCP's and 
IMF's with Z < 20 in « 40% of the forward solid angle 
(plus a reduced sampling in the backward hemisphere). 
The phoswiches were two-element modules (fast plastic 
BC404 + CsI(Tl)) or three-element ones (BC404 + slow 
plastic BC444 + CsI(Tl)), coupled to fast phototubes. 
The fast scintillators had been carefully machined down 
to 200 /im in our workshop, with thickness uniformity 
better than 5% as required for Z identification up to Z w 
20; the thickness of the BC444 was 5 mm and that of 
the CsI(Tl) 30 or 50 mm. In the three-element modules, 
the lack of dead layers and the good optical transmission 
(obtained by coupling the two plastic elements with the 
heat-pressing technique) allowed to identify IMF's with 
low thresholds («3-10 ^MeV for Z=2-20). The measure- 
ment of the time of flight allowed to directly obtain the 
velocities of LCP's and IMF's, without time consuming 
and tricky energy calibrations of the scintillators. A de- 
tailed description of the apparatus can be found in |lj| . 

The data presented in this Letter are focused on binary 
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events, where only two major fragments are detected in 
the PPAD's while the coincident LCP's and IMF's hit the 
phoswich detectors. Using the Total Kinetic Energy Loss 
(TKEL) as an ordering variable, it is possible to select 
events with increasing impact parameter up to grazing 
collisions. At variance with respect to all available Att 
detectors, our low-threshold apparatus makes it possible 
to measure both the PLF and the TLF also for events 
with TKEL< 400 MeV. Reconstructed primary quanti- 
ties were obtained from the KCM with 2-body kinemat- 
ics, not only for LCP, but also for IMF emission |l6[ |. 

The left panels of Fig. |l| present, for TKEL=240-400 
MeV, the experimental yield of p, d, He and IMF (Z=3- 
7) in the V perp - V par plane, V par and V perp being the 
center-of-mass velocity components parallel and perpen- 
dicular, respectively, to the separation axis of the two ma- 
jor fragments. The dotted lines show the velocity thresh- 
olds due to the thickness of the thin plastic scintillator. 
The data have been corrected for the finite geometry of 
the apparatus because of its large acceptance and 
axial symmetry, the correction is largely independent of 
the emission pattern. In the backward lab-hemisphere 
{Vpar^ —40 mra/ns), owing to the reduced detector cov- 
erage, the correction is not as effective as forwards. 
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FIG. 1. Left: experimental yields in the V per p-V pa r plane 
for p, d, He and IMF's (Z = 3 - 7) in the system 116 Sn + 
93 Nb at 29.5 AMeV, for TKEL =240-400 MeV (corrected for 
the setup geometry). Level spacing is logarithmic, dashed 
lines indicate velocity thresholds. Right: same results for the 
simulated statistical emission from the hot reaction partners. 

In order to better put into evidence the major fea- 
tures of the experimental data, the right panels of Fig. |l| 
show the corresponding yields obtained with the simu- 
lation of a pure statistical emission from fully acceler- 
ated fragments, filtered with the setup acceptance and 
then corrected as the experimental data. The excitation 
energy of the fragments was obtained from TKEL with 



an "equal energy" sharing, according to ]lq| , while first 
guesses of the parameters of the evaporation step were 
deduced from calculations with the code Gemini |l7j . 

The experimental emission pattern for protons does 
not differ very much from that expected for a sequen- 
tial evaporation, but with increasing particle mass the 
mid-velocity yield becomes increasingly important until 
it actually exhausts most of the IMF intensity. 

The mid-velocity particle multiplicities were obtained 
with a procedure similar to that outlined in [||. The 
simulation was tuned so as to reproduce the experimen- 
tal data in the velocity region corresponding to forward 
emission from the PLF (10° < 9 PLF < 40° in the PLF 
reference frame). Assuming that the simulation prop- 
erly mimics the whole statistical emission for PLF and 
TLF, the mid-velocity yield is obtained as the difference 
between the total experimental emission and the corre- 
sponding estimation of the statistical component. 
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FIG. 2. Experimental multiplicities of p, d, t, He and IMF's 
(Z=3-7) against TKEL in 116 Sn + 93 Nb at 29.5 AMeV. Full 
squares (circles) refer to the mid- velocity (evaporative) com- 
ponent; open symbols are for the system Xe + Sn at 32 AMeV 
[0. Lines are to guide the eye. On top, correspondence be- 
tween TKEL and impact parameter (or centrality binning 
of 0), estimated from the QMD code CHIMERA Q. 

The deduced multiplicities are shown as a function of 
TKEL in Fig. |, separately for p, d, t, He and IMF. Full 
circles refer to the statistical evaporation, with uncer- 
tainties (due to detector thresholds and determination of 
the TLF emission) usually smaller than the symbol size. 
Full squares show the multiplicity of mid-velocity parti- 
cles. The larger the mid-velocity component, the more 
reliable its extraction is; the largest errors are for protons, 
owing to the presence of a large statistical component. 
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It is worth noting that, in the mid- velocity emission 
of peripheral events, the mass removed by the IMF's (as- 
suming A w 2Z) is comparable to that removed by LCP's 
(while it is much smaller in case of evaporation) . More- 
over, as already qualitatively seen in Fig. ^ the mid- 
velocity component of IMF's greatly overcomes the sta- 
tistical emission, by more than a factor of 20 for the most 
peripheral collisions. To our knowledge, this extremely 
enhanced emission of mid-velocity IMF's has never been 
directly observed before. In these events the nuclear mat- 
ter tends to break apart in intermediate mass fragments, 
this process successfully competing with LCP emission. 

The comparison of the multiplicities presented in Fig. |^ 
with those of Xe+Sn (open symbols) at similar bom- 
barding energies 0, requires to find a correspondence 
between our ordering variable (TKEL) and that (trans- 
verse energy of LCP's, E trans i 2 ) used by the INDRA col- 
laboration to estimate the impact parameter. Following 
ref. §, we used the same QMD code CHIMERA flf] as 
an event generator for our reaction 116 Sn + 93 Nb. In 
simulated events, analyzed as the experimental ones, the 
reconstructed TKEL is narrowly correlated with the im- 
pact parameter ||o) ; details are given in |l5| . Here, on top 
of Fig. |[ we synthetically draw the obtained impact pa- 
rameter scale, the notation underneath (BIN1, BIN4) 
corresponding to the binning of 0. It has to be noted 
that the range of impact parameters probed by our ex- 
periment extends to significantly larger values than those 
accessible with INDRA, the results of both experiments 
being in good agreement in the common region. 

Addressing the question whether the production mech- 
anism of mid- velocity IMF's is mainly of dynamical or 
statistical nature is beyond the scope of this letter. How- 
ever, the peculiar emission pattern observed in the very 
peripheral collisions (where it is almost free from signif- 
icant contaminations due to evaporation from PLF and 
TLF) gives information on the time-space configuration 
of the emitting system. Namely, it is possible to infer the 
time scale and the geometrical configuration at the end of 
the reaction from Coulomb trajectory calculations repro- 
ducing the observed emission pattern. Although qualita- 
tive descriptions of the IMF emission based on Coulomb 
trajectory calculations are found in the literature (see, 
e.g. pi]]), quantitative conclusions are not frequent. 

A neck-region emission was assumed, with the IMF's 
emitted from a "neck" of excited nuclear matter con- 
necting the two interacting nuclei. Schematically, the 
Coulomb trajectory calculation was implemented in our 
Monte Carlo simulation as follows. The two main frag- 
ments were assumed to be in an initial configuration (cor- 
responding to the selected peripheral impact parameter) 
at a distance given by a tuning parameter (d sop ) of the 
calculation. The to-be-emitted IMF, sampled from a re- 
alistic distribution, was located in between the two major 
fragments, with a random "thermal-like" kinetic energy 
described by another tuning parameter (E n ). The ini- 



tial separation velocity of the two main fragments was 
that appropriate for producing the correct asymptotic 
TKE. The equations of motion were integrated starting 
from this simple initial configuration in the three-body 
phase space, then the experimental filter was applied and 
the V r 



perp 



V P ar plots were built. Fig. ||a presents the 



obtained emission pattern for impact parameters corre- 
sponding to the experimental data in the bin TKEL= 
240-400 MeV. The main effect of the two close-lying frag- 
ments is to focus the emitted IMF in the transverse di- 
rection, with the shape of the V perp distribution mainly 
depending on the assumed thermal-like energy. In fact, 
for peripheral collisions, the heavy fragments fly swiftly 
apart, with no significant energy transfer to the IMF's 
if they are emitted nearly at rest in the CM system. 
From the comparison with the experimental correlation 
(Fig. |^d), one sees that the simulated mechanism of neck- 
region emission significantly contributes to the experi- 
mentally observed mid- velocity emission, in particular in 
the high transverse momentum region, but there is no 
way to reproduce the two Coulomb-like wings which in 
the experimental data start at mid-rapidity and rapidly 
fade away when moving forward (backward) in the PLF 
(TLF) emission frame. 
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FIG. 3. V perp -Vp a r plots for IMF's (Z=3-7) at TKEL= 
240-400 MeV: (a) Coulomb trajectory calculation for 
neck-region emission with (E n ) =25 MeV and (b) for sur- 
face emission with (E 3 ) =10 MeV, (c) superposition of the 
two contributions and (d) experimental correlation. Projec- 
tion of experimental (full histogram) and simulated (dotted 
histogram) yields on (e) the V par and (f) the V pe r P axis. 

It was therefore assumed that a second contribution 
exists, similar to the so-called fast oriented fission [jll 22 
where the very asymmetric fission-like decays of A sa 
100 nuclei proceed through nearly aligned configurations 
(light fission fragment in between the heavier one and the 
other non-fissioning reaction partner) with characteristic 
times of the order of 10 -21 s « 300 fm/c. We performed 
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Coulomb trajectory calculations assuming a fast surface 
emission of IMF's from the contact regions on the surface 
of the two flying-apart primary fragments. The IMF was 
emitted, again with a "thermal-like" energy (E s ), from 
a point of the surface randomly spread within a few fm 
from the separation axis. The emission pattern obtained 
for such emissions is presented in Fig. ^b. This surface 
mechanism reproduces the wings of the experimental dis- 
tribution (Fig. ||d), but not the large transverse energies: 
in this case the distribution is mainly determined by the 
Coulomb field of the emitting fragment. 

In order to satisfactorily reproduce the experimental 
data, both contributions must be included. The best 
agreement with the experiment has been obtained for 
a ratio between surface and neck-region emission of 0.7 
and with values of (E s }= 10 MeV and (E s )= 10 MeV: 
the result is shown in Fig. ||c. The projections of the ex- 
perimental data and the simulation (continuous and dot- 
ted histograms, respectively) on the V par (Fig. ||e) and 
V perp (Fig. ||f) axes allow to better perceive the good 
quality of the obtained agreement. In both projections, 
the simulation presents a sizeable sensitivity to the as- 
sumed random energies (E n ) and (E s ). In fact, a change 
of 5 MeV in these parameters worsens in a significant 
way the agreement. An agreement of the same quality 
is found for more peripheral events (TKEL<240 MeV), 
while for more central TKEL-bins the comparison be- 
comes increasingly blurred by the growing contribution 
of the statistical emission. The overall trend is that both 
the random energy (E n ) of the neck-region emission and 
the spread of the emission point in the surface emission 
increase with increasing centrality. 

For peripheral events, the experimental IMF emission 
appears to be compatible with the formation of a neck- 
like structure which then decays by a prompt emission 
from the neck region itself or by a successive emission 
from the surfaces of the separating nuclei. For example, 
the neck-region emission may be due to a multiple neck 
rupture, while the surface emission is suggestive of a sin- 
gle neck rupture p3f , leaving one of the main fragments in 
a rather deformed shape, in the vicinity (or even beyond) 
its saddle shape. In the simulation these two mechanisms 
differ in terms of the initial configuration: the neck-region 
emission contributes at small distances d sep compatible 
with the presence of a connecting neck, the surface emis- 
sion at larger distances. According to models quanti- 
tatively describing the formation of the neck in nuclear 
collisions or fission p^j , the maximum length of the neck 
-discriminating between the two mechanisms- was esti- 
mated to be 12 fm. From the relative kinetic energy at 
separation and from the envisaged geometrical configu- 
ration, it is possible to estimate the time scale: a value of 
the order of 60 fm/c is obtained as a discrimination be- 
tween neck-region and surface emission. Moreover, from 
the observed anisotropy and using realistic values of the 
angular momenta and moments of inertia of the nuclei, 



one can exclude for the surface emission times signifi- 
cantly larger than a few hundreds of fm/c. 

In conclusion, peripheral collisions are characterized by 
an important emission of IMF's at mid- velocity, success- 
fully competing with LCP's. Indications for a fast energy 
dissipation process of locally excited nuclear matter may 
be inferred from the estimated thermal-like energies nec- 
essary to reproduce the data. The emission pattern of 
IMF's is compatible with the coexistence of two mecha- 
nisms of prompt or fast emission: the first one related to 
the formation and prompt break-up of a neck of highly 
excited nuclear matter and the second one, at a somewhat 
later stage, characterized by a localized emission from the 
possibly highly deformed flying apart fragments. 
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